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I. INTRODUCTION 

About one and half decades ago, spin models defined 
on hierarchical lattices received much attention in the lit- 
erature US S S S B S II in • In general, to construct 
a hierarchical lattice we first start with an unit, which 
may be a bond or a cell, and then proceed a given type 
of bond- or cell-decoration iteratively to the infinite limit. 
Thus, a hierarchical lattice has fractal structure, and the 
thermodynamic limit for a physical system defined on 
a hierarchical lattice is well defined. Hierarchical spin 
models attract researchers' interest mainly due to two 
reasons; Firstly theses models are exactly solvable in the 
context of the Migdal-Kadanoff renormalization scheme 
[13, [13. Secondly owing to the inhomogeneity in the co- 
ordination number of lattice sites some particular prop- 
erties revealed from the models may provide insights to 
inhomogeneous systems such as random magnets, poly- 
mers, and percolation clusters |l2j| . 

Q-state Potts model defined on a diamond hierarchical 
lattice is one example. Starting with a bond, a diamond 
hierarchical lattice is obtained by replacing bonds by di- 
amonds iteratively to the infinite limit. There exists a re- 
markable richness of phenomena for the model in the ab- 
sence of external fields. In particular, the limiting set of 
the partition function zeros in the complex temperature 
plane, also referred as the Fisher zeros, are essentially 
the Julia sets associated with the rational map defined 
by renormalization transformation Q , and the Julia 
set possesses multifractal structure for g > [l^ 113. Il5l| . 

The interest about the loci of partition function zeroes 
has been raised after the classical works of Yang and Lee 
on regular lattices 0, 0|. After the remarkable Lee- 
Yang circle theorem. Fisher studied partition function 
zeros in the complex temperature plane and showed that 
the distribution is an unit circle in the sinh (2 J/fesT) 
complex plane for the two-dimensional zero-field Ising 
model on simple square lattice jl^. Since then, the dis- 
tributions of Fisher zeros of the Ising model with isotropic 
or anisotropic couplings on a variety of classic planar lat- 
tices have been investigated [T9II20II21I . Recently, Lu 
and Wu completed the picture by calculating the density 
of zeros for two-dimensional Ising model in zero field as 
well as in a pure ir nag inary field i7r/2 on a variety of clas- 
sic planar lattices [23j. In principle, by knowing the zeros 
of the partition function and the corresponding densities, 
we may deduce all the thermodynamic characteristics of 
a system. For example, the density of the zeros near the 
critical point can be used to extract the critical expo- 
nents ,20, JlJ , and the logarithmic singularity of the spe- 
cific heat for the two-dimensional zero-field Ising model 
is the result of the linearly vanishing density of the zeros 
near the real axis [isL |23| . 

In this paper we study the distributions and the den- 
sities of the Fisher zeros of the zero-field Ising model on 
square lattices with diamond-type bond-decorations, re- 
ferred as diamond-decorated Ising model (DDIM). The 
lattices used in the model are constructed by starting 



with a simple square lattice, and then by implemend- 
ing diamond-type bond-decorations to each bond itera- 
tively to any desired degree. For DDIM, there exists a 
well-defined thermodynamic limit for any finite degree of 
decorations, and each primary bond becomes a diamond- 
hierachical lattice used in the diamond-hierarchical spin 
model for the limit of infinite decorations. In our pre- 
vious work, the properties of ferromagnetic phase tran- 
sitions of DDIM have been investigated extensively for 
finite as well as infinite decoration-levels Here we 

concentrate on the distribution and the density of the 
Fisher zeros. 



Similar analyses on the distribution of the Fisher ze- 
ros have been carried out for triangular type Ising lat- 
tices with cell-decorations [2^ . These lattices possess the 
Sierpihski gasket as the inherent structure for a primary 
triangle in the limit of infinite decoration-level. The re- 
sults indicate that the zeros distribution for the infinite 
decorated lattices coincides with those for the model de- 
fined on the Sierpihski gasket, and the distribution of 
zeros appears to be an union of infinite scattered points 
and a Julia set called the Jordan curve, and the scattered 
points are bounded by the Jordan curve. Note that the 
Jordan curve is a quasi-one-dimcnsional circle with the 
Hausdorff dimension equal to one. It is also well known 
that the limiting set of the distribution of the zeros of 
the diamond-hierarchical Ising model (DHIM) is a Julia 
set which owns a multi- fractal structure 0,0,01^3.^3- 
The Julia set, which is a bounded planar distribution 
with the Hausdorff dimension greater than one, is quite 
different from the Jordan curve. Based on these obser- 
vations, we may expect that the Julia set occurring in 
DHIM may also appear in DDIM, and thence we may un- 
derstand the formation of the multi-fractal structure in 
the Julia set by studying the variation of the distribution 
and the density of the zeros in the course of increasing 
the decoration- level of DDIM to the infinite limit. 



This paper is organized as follows: In Section II, we 
briefly describe how to deduce the exact expression of free 
energy via the construction of the exact renormalization 
map of temperature between two successive decoration- 
levels and the use of the known results of the Ising model 
on simple square lattice. In Section HI, we study the 
distribution of the Fisher zeros and exhibit the change of 
the distribution pattern as the decoration-level increases. 
In Section IV, we determine the density of the zeros for 
the first two decoration-levels by using the results for 
the case of simple square lattice, and the densities for 
higher decoration-levels are given qualitatively by con- 
jecture. In Section V, we discuss how the Jufia set arise 
in the limit of infinite decoration-level and characterize 
its global mutifractal structure. Finally, Section VI is 
preserved for summary and discussion. 
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II. FREE ENERGY 

We construct the exact expression of the free energy 
of DDIM with an arbitrary decoration-level n in this sec- 
tion. A simple square lattice with diamond-type bond- 
decorations up to the level n is referred as an n-lattice. 
Then, simple square lattice itself is 0-lattice, and its con- 
necting bonds between any two nearest neighbors are 
named as 0-bonds. We denote the total site-number and 
bond- number of 0-lattice as ris and rib with rib — 2ns. 
An 1-bond is formed by replacing a 0-bond by a diamond 
which consists of four 0-bonds. Starting with a 0-bond, 
after the n-fold iterative replacements of 0-bonds with 
1-bonds, we obtain an n-bond which has site-number 
gin) ^ 2(4" + 2)/3 and 0-bond-number = 4". 

An n-lattice is formed by replacing all 0-bonds of a 0- 
lattice with n-bonds. For an n-lattice, the average site 
and bond numbers per primary square of 0-lattice are 



TVi"^ = 25'(") - 3 and NI""' = 2B(") respectively. The 
construction procedure is schematized in Fig. 1. 

The general form of the partition function for DDIM 
defined on an n-lattice reads 



E 



Y[ exp (rjaiaj) 



(1) 



where the sum is over all bond-connected pairs 
of n-lattice, and the Ising spin takes two possible val- 
ues (Ti = ±1. Here we consider uniform ferromagnetic 
couplings characterized by the coupling strength J, and 
the dimensionless coupling parameter rj is defined as 
r] ^ J/ksT. 

To calculate the partition function of Eq. (1) for 
an arbitrary n-lattice, we use the bond-renormalization 
scheme in evaluating the Boltzmann factors associated 
with an n-bond. The details of the derivations are given 
in Ref. p^ . and we briefly summarize the results in the 
foUowings. 

The Boltzmann factor associated with the 0-bond, de- 
noted by sj'^-'^,, is given as exp (rja^a^), and it can be 



written as 



B 



(0) 



cosh (ry) + ct^CTi, sinh {r]) . 



(2) 



There are decorated (S'(") - 2) -sites for an n-bond. The 
Ising spins defined on the decorated sites couple only to 
those belonging to the same n-bond, and we refer them 
as inner spins. Then we may define the Boltzmann fac- 
tors associated with an n-bond, B 



as the result of 
taking the sum over the inner spins for the product of all 
Boltzmann factors associated with the n-bond. 



B 



(n) ^ 



(s(")-2) 



Here the two subscripts, and v, denote the two pri- 
mary sites before decorations, the front factor is added 
for the normalization of the sum, and the sum is over the 
_ 2)-inner spins. By substituting the expression of 
Eq. (2) into each Boltzmann factor of Eq. (3), we have 



= ^'"^ (t^^")) + a^cj, sinh (t?'"' 

for n > 1, where the function i?^"^ (77) is given as 



n 



fe=i 



and 77'-'^'' is determined by the recursion relation. 



exp 1 77 



cosh 



(4) 
(5) 

(6) 



with the initial condition 77*^*^^ — 77 for 1 < fc < n. 

It is well known that the corresponding free energy per 
bond per ksT oi Eq. (1) for the case of 0-lattice can be 
written as 



J(0) = _ 



- In sinh (277) 

/ 2^ X ^ 

-9(0,0)], 



sinh (277) 



sinh (2?7) 



with 
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cos W -h COS ( 



(7) 



(8) 



By observing that up to a factor i?*^"-' (77) the effective 
Boltzmann factor of an n-bond posesses the same form 
as that of a 0-bond, W6 can express the free energy density 
of an n-lattice as |2j 



/ 



(n) _ An) _ 

— Jl 



1 



4B(") 

r2v 



In sinh ( 2ri 



^(") 



de 

1 



In 



sinh 



(277(")) 



sinh (2r/( 



n)) 
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(9) 



where /jj'-* is the contribution from the factor i?*^"^ (77), 



/l.'' = -ilni?(")(77). 



which can be expressed as 



fe=i 



(10) 



(11) 



Ea^-^\^f^ ^ ^ ^ _L^^A'i> by using the recursion relation of Eq. (6). Note that the 

exp [77 (cr^O-a + C^aCTb + ... + 0-sCry(>]l) („) . ^ w 



f'^^ part exists only for n > 1. 
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III. FISHER ZEROS 

The partition function zeros of a 0-lattice in the ther- 
modynamic limit can be obtained by setting the argu- 
ment of the logarithm in the free energy density of Eq. 
(7) equal to zero |2^. It is known that the zeros may lie 
on the unit circle, |sinh(2ry)| = 1 or on two circles, 
|tanh(ry) ± 1| = y/2 0, depending on the variable used 
for the complex temperature plane. In this paper, we 
study the distribution and density of the Fisher zeros in 
the complex tanh (rj) plane. The basic features appearing 
in the complex sinh (2r]) plane are essentially the same as 
those we obtain in the complex tanh (77) plane. 

By observing the free energy density of Eq. (9), in the 
thermodynamic limit we can obtain the zeros distribution 
of an n-lattice from the solutions of two conditions. 



sinh 



(2,(")) 



sinh {2ri<''' 



-ei 



= 0, 



and 



cosh 



(2,7(^-1; 



for fc = 1,2, 



(12) 



(13) 



Note that the latter can also be viewed as the condition 
of the zeros for the Ising system defined on an n-bond, 
and such a system becomes DHIM in the limit of infinite 
n. 

^From the result of two circles for the partition func- 
tion zeros of a 0-lattice in the complex tanh (77) plane, we 
may conclude that the solution of Eq. (12) is two circles 
in the complex tanh (r;'-"-') plane. 



tanh(77("M ± 1 



= x/2. 



(14) 



The two circles intersect at two points, i and —i. As 
shown in Fig. 2a, due to the intersections there is a ring 
contained in the two circles. For the purpose of identifi- 
cation, we refer the circles as two n-cycles and the ring as 
n-ring. Note that in displaying the distributions of zeros, 
we always use bold-faced curves for the right n-cycle and 
its descendants, and gray curves are for those from the 
left 7^-cycle. This distribution also appears to have the 
symmetry of inversion about the center tanh (77^"^) = 0. 

For the purpose of comparing the zeros distributions 
among different n-lattices, we have to bring the zeros to 
the complex plane of an unique variable chosen to be 
tanh (77). To achieve this, we notice that the recursion 
relation of Eq. (6) can be rewritten as 



tanh (77*") 



2 [tanh (77("~i))]^ 
1 + [tanh (7/("-i))]' 
which has the inverse map given as 



(15) 



tanh (77^""^^ 



1 ± 



± 



^Jl- [tanh (7/(") 



1/2 



tanh (7/(")) 



Then, starting with the two 7i-cycles in the complex 
tanh (77^"^) plane, we can obtain the corresponding zeros 
distribution in the complex tanh (77) plane by performing 
the Ti-fold backward iterations provided by Eq. (16). 

After the first backward iteration, we show the resul- 
tant distribution in the complex tanh77("^^^ plane in Fig. 
3a where the points indicated by crosses are the preini- 
ages of the map of Eq. (15) for the centers of two 77- 
cycles, 1 and —1. The results indicate that each of the 
ri-cycles shown in Fig. 2a splits to two closed curves re- 
ferred as (77 — l)-cycles. There are 8 intersection points 
between the descendants of the right 7i-cycle and those 
from the left n-cycle, and the loci of the intersections are 
determined by the inverse maps of the points, i and —i, 
which are the intersections of two 7T.-cycles. There are 
four rings, referred as (77 — l)-rings, caused by the inter- 
sections. Note that the distribution shown in Fig. 3a are 
just the zeros distribution in the complex tanh (77) plane 
subject to the condition of Eq. (12) for an 1-lattice. 

Proceeding with the inverse map given by Eq. (16) 
from the complex variable tanh (77^""^') to tanh (t;^""^^) 
for the four {n — l)-cycles, we obtain 12 (77 — 2)-cycles 
as shown in Fig. 4a. The {n — 2)-cycles contain 16 
{n — 2)-rings caused by the 32 intersections among the 
cycles. The intersections are again the preimages of the 
map of Eq.(15) for the loci of the intersections among the 
{n — l)-cycles. 

Continuing with this procedure, we show the distri- 
bution of zeros in Fig. 5 for n = A and Fig. 6 for 
77, = 8. In general, the zero distribution for a 77-lattice 
in the complex tanh 7/ plane, subject to the condition 
of Eq. (12), is the union of [2 + 2 (4" - 1) /3] 0-cycles 
which have 4"+^/2 intersections, and these intersections 
yield 4" 0-rings contained in the 0-cycles. The 0-cycles 
can be divided into two classes: one consists of the the 
descendants of the left 7X-cycle with [2-^2 (4"-^ - l) /3] 
members, and the other has 4"/2 members which are the 
descendants of the right 77-cycle. The intersections only 
occur between two 0-cycles belonging to different class. 
The distribution always maintain the inversion symmetry 
about the center tanh 77 = 0. 

For the condition of Eq. (13), it can be rewritten in 
terms of the tanh77'^"^ variable as 



tanh 77 



(fc-i) 



±i for fc = 1, 2, 



(17) 



(16) 



Then, for the case of 77 = 1 the two zeros, i and —7, in the 
complex tanh 77 plane are just the intersection points of 
two 77,-cycles in the complex tanh(77'^^^) plane. Proceed- 
ing to the case of 77 = 2, we obtain, besides of the previous 
two zeros, 8 more points in the complex tanh (rj) plane 
from fc = 2 in Eq. (17). These additional points are the 
preimages of the renormalization map for the original two 
points, and they are the intersection points among four 
{n — l)-cycles in the complex tanh (77'^""^^) plane. By in- 
duction, we may conclude that the zeros in the complex 
tanh (77) plane obtained from the condition of Eq. (17) 
are 2 (4" — 1) /3 scattered points which are the union of 
the intersection points among fc-cycles in the complex 



4 



tanh(77^'^^) plane for 1 < k < n. in the complex tanh?] 
plane for an n-lattice. 

Thus, we can describe the distribution pattern of the 
Fisher zeros of DDIM with the decoration-level n in the 
complex tanhrj plane as follows. There are 2 (4" — 1) /3 
scattered points given by Eq. (17). In addition, there 
are [2 + 2 (4" - 1) /3] 0-cycles with 2 • 4" intersections 
obtained from Eq. (12). 

Among all the zeros we obtained in the above, as a 
consequence of the Lee- Yang theorem 0, 0| , the bulk 
critical points correspond to the zeros falling on the phys- 
ical region. The physical region of the variable tanh (rj) 
is < tanh (77) < 1 for ferromagnetic couphngs 77 > 0. 
This implies that the variable tanh (77^"-*) also takes the 
range < tanh (77*^"^) < 1 as the physical region for any 
n value. 

/^From the solutions of the conditions of Eqs. (12) and 
(13), we know that there is only one zero in the physical 
region. This zero belongs to the right 7z-cycle and locates 
at 



tanh 



(18) 



with h^^^f — \/2 — 1 for arbitrary decoration-level n. 
Here, for convenience, we use the notation, =, to de- 
note the equality established only at the phase transition 
point. Note that the /i^"^ value is just the reduced critical 
temperature, tanh ( J/fc^Tc), of the ferromagnetic phase 
transition for the square Ising model, i.e. tt, = 0. 

To find the locus of the zero specified by Eq. (18) in the 
complex tanh 77 plane, we can continuously use the inverse 
map of Eq. (16) to obtain the equivalent expression of 
Eq. (18) as 



tanh 



with 



1/2 



(19) 



(20) 



for 1 < fc < 77,. Note that in obtaining Eq. (20) for the 
critical value of tanh 7;*^""'^') we have used the constraint 
< tanh77("^''') < 1. 

Thus, the zero, , in the complex tanh 77^"^ plane cor- 
responds to the zero, in the complex tanh 77 plane, 
and the /i^") value is just the reduced critical tempera- 
ture tanh ( J/fc^Tc) of the ferromagnetic phase transition 
for DDIM with the decoration-level n. The sequence of 
/i*-"^ decreases as n increases, and the ft,^"^ value in the 
limit of infinite n is given by the asymptotic value of 
the sequence of /i'"-'. For the recursion relation of Eq. 
(15), there are three fixed points, one repellor locating 
at 0.5437..., and two attracors at and 1. Since the ft.'^"-' 
value is obtained from h^'^^ via the 71-fold backward iter- 
ations given by Eq. (20) and the attractors (repellors) of 
the map become the repellors (attractors) of the inverse 



map, we may conclude that the /i^"-* value is the locus of 
the repellor of Eq. (15), /i^") = 0.5437..., for the case of 
infinite 77. 



IV. DENSITY OF ZEROS 

The density of two 7i-cycles in the complex tanh77'^") 
plane has been determined by Lu and Wu Based 
on this result, we determine the density of the zeros of 
DDIM in this section by performing proper transforma- 
tions. First, we describe the result of Lu and Wu briefly 
in the following. 

The two cycles of Eq. (14) are written as 



tanh77(") ± 1 



r exp fia^") 



(21) 



where r (o;'^"-') = \/2 is the radial distance from the cen- 
ter coordinate. By considering an M x 2N simple-quartic 
lattice with Brascamp-Kunz boundary condition, we in- 
troduce the zero density, g± (a'-"-'), which satisfies the 
normalization condition 



2tx 



(22) 



such that the number of zeros in the interval 
[q,(»)^q,(") +dck(»)] is 2MiV5± («(")) da(") for the left 
(-f ) and right (— ) n-cycle respectively. The zero density 
is given as 



("'•")-»-('-"'■")- (3 



= .9 

xK{x) 



1 - V2cosa(") 



2V2cosa(") 
(23) 



where 



2|sina(")| (\/2 - cos a^")) 



2V2 



cos a 



(") 



(24) 



and K (x) is the complete elliptic integral of the first 
kind, 



K{x) 



dt- 



1 



\/\ ~ sin'^ t 



(25) 



The density of g_(-(a''"'), which has the symmetry 
5+ = g+ (27r - a'")), is plotted in Fig. 2b for the 

range < a*^"^ < tt. Here we also specify the radial dis- 
tance of a zero, r («(")), in the right vertical scale. For 
small q;'-"\ the result of Eq. (23) has the linear behavior 
as 



'3±2V2 



v(") 



(26) 



Note that the the density g+ (a^"-*) of the zeros near the 
the point a'"' = is the same as the density g- (a^"-*) 
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near the point which is the ferromagnetic phase 

transition point of the bulk system, and this hnearly van- 
ishing density of the zeros near the bulk transition point 
leads to the logarithmic singularity of the specifi heat. 

To find the corresponding density in the complex 
tanh (?7("~^)) plane, we first write the (n — l)-cycles as 



tanhr^ 



(n-l) 



+ r (a("-i))exp(ia("-i)). (27) 



Here the center coordinates, Zq, are chooscn to be the 
preimages of the map of Eq. (15) for the center co- 
ordinates of n-cycles, tanhTy^"^ = 1 and —1, and this 
leads to zq = 1, —1, i, and —i for the center coordi- 
nates of 4 (ri — l)-cycles. The zero density of the (n — 1)- 
cycle, specified by the center coordinate zq, is denoted as 
Qzo (q;("~^^), and we have the relation, 



51 



(a(-^) - I) , (28) 



= 9- 



for the distribution shown in Fig. 3a. Thence, we deter- 
mine the density gi (a*^"~^^), and the densities of other 
(n — l)-cycles are followed from gi (a^""^^) according to 
the above relation. 

Because that the inverse map given by Eq. (16) is one 
to four and the (n — l)-cycle of = 1 is a descendant of 
the left n-cycle, we can express the density gi (a^"~^^) 
as 
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(a("-i)) = 



9+ (a 



in)) 



(29) 



To determine the transformation Jacobian, 
|rfa^"^/da^"~-^^|, we first notice that the (ra — l)-cycles 
of = ±1 are the solutions of the equation, 



tanhT?^"-! 
-2V2 tanh??^"-^) 



(tanhry^"-!)) V (tanhry^""!)*)^ 



+ 1 



= 0, 



(30) 



where tanh the complex conjugate of 

tanh?]^""'^^. This result is obtained by substituting Eq. 
(15) into Eq. (14) for the left n-cycle. By further substi- 
tuting Eq. (27) with zq = 1 into Eq. (30), we obtain 



4 cos 



— r 
0. 



4V2 



cos 



(4 - 2V2) 



2V2 



(31) 



This equation can be solved numerically to obtain 

r (a*^""^)) and dr/da'^^^^^h Moreover, the relation be- 
tween and a'""^) has been specified by the map of 



Eq. (15). We substitute Eqs. (21) and (27) into Eq. (15) 

to obtain 



l + \^exp(ia(")) 



2[l + r (a("-i)) exp(m("-i))] " 



4 ■ 



l+[l + r («("-!)) exp(ia("-i))] 

(32) 

By differentiating this equation with respect to 
q;("-i) and by using the known values of r(a*^"~^-') 
and dr/da^"~^\ we can determine the derivative, 
da^"") / da'^"~^\ and then to obtain the density 
gi (a("-i)) according to Eq. (29). 

The numerical result of the density g^ (a^""^^) is 
shown in Fig. 3b for the range < < TT with 

the radial distance of a zero, r (a^""^)), specified in the 
right vertical scale. Our results indicate that when the 
complex plane changes from tanhry^"^ to tanh 77'^"^"'^', the 
distribution density oscillates more rapidly with the peak 
number increasing from 2 to 4 for half cycle. The locus of 
the zero corresponding to the ferromagnetic phase tran- 
sition point moves from a*^"^ = of the left n-cycle to 



(n-l) 



a 

near to 
as 



TT of the (n — l)-cycle of zq = 1. For the zeros 
^) = TT, the density has the linear behavior 
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rn 









with 



^-4 



0.1529. 



(33) 



(34) 



v(n-l) = 



This linear behavior, again, gives the logarithmic singu- 
larity of the specific heat. 

To extend the calculation of density to (n — 2)-cycles, 
we may divide the 12 (n — 2)-cycles shown in Fig. 4a into 
three classes named as class /, //, and ///, according to 
the decreasing order from the longest to the smallest in 
the magnitude of the circumference of the circles. Then, 
there are four members in each class, and all members of 
class I are the descendants of the left n-cycle and those 
beonging to class // and /// are from the right n-cycle. 
Similar to the case of (n — l)-cycles, we can write 

tanhry("-2) =zi+r (a^"-^)) exp(ia("-2)) (35) 

for the zeros of (n — 2)-cycles, and the center coordinates, 
zi , are chosen to be the preimages of the map of Eq. (15) 
for the center coordinates of (n — 2)-cycles, 1, —1, i, and 
—i. Then, we can express the densities of the (n — 2)- 
cycles as 



9k 



91 (a 



("-1)1 



and 



9iy'' 



(a 



(n-l)) 



da("-2) 



(36) 



(37) 
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respectively, where the superscript, /, //, or ///, is used 
to denote the class to which a. {n — 2)-cycle belongs, and 
the subscript, zi, is used to specify a (n — 2)-cycle in 
the given class. Since the members belonging to the 
same class are the same up to a global rotation, we only 
need to determine the zero density of a (n — 2)-cycle for 
each class. The cycles of zi — 1, \/l + \/2e*'^/^. 



and 



-l + -\/2e™/^, belonging to class /, //, and /// re- 
spectively, are choosen for the calculation of the density 
of the respective class. 

The numerical method of calculating the correspon- 
dence between o;'-"^^' and a*^"^^' and the Jacobians, 
|(^Q,("-i)y'(^Q,("-2) I ^ are exactly the same as we did in the 

last case. The results of r (a*^""^-*) (right vertical scale) 
and g (a'-""^-') (left vertical scale) are shown in Figs. 4b, 
4c, and 4d respectively for the three cycles. These results 
indicate that the rapidity of oscillation in the distribu- 
tion density of the cycles of class / increases as the peak 
number doubles with respect to the last case, while the 
peak number remains to be the same for the cycles of 
class // and ///. The zero corresponding to the criti- 
cal point of ferromagnetic phase transition moves from 
Q,("-i) = TT of the (n — l)-cycle of zq = 1 to a^""^) = tt 
of the (n — 2)-cycle of zi = 1 of class /. The density of 
the zeros near to this locus has the linear behavior. 



(7r + a("-^))=<52 



3 + 2V2 



v("-2) 



(38) 



with 



So = - 



da'^' 



0.0365. 



(39) 



^From the densities obtained in the above, we may 
conjecture qualitatively the density of the 0-cycles with 
n = 4 shown in Fig. 5 as the following: There are 128 
members belonging to the descendants of the (n — 2)- 
cycles of class // and /// as displayed by bold-faced 
curves in Fig. 5. The 16 members of the 128, appearing 
in the outermost of Fig. 5, are similar to a (n — 2)-cycle 
of class //, and the rest have a similar shape as a (n — 2)- 
cycle of class ///. Up to an overall reduction factor, the 
densities of the 16 members have the same oscillation 
pattern as that shown in Fig. 4c and the densities for 
the rest members possess the same oscillation pattern as 
the one shown in Fig. 4d. For the 44 descendants of the 
(n — 2)-cycles of class / displayed by gray curves in Fig. 
5, the corresponding density reduces but oscillates more 
rapidly with 32 peaks in half cycle in comparing with the 
one shown in Fig. 4b. 

Continuing with this procedure in finding the density, 
we obtain 2"+^ peaks in a tt period for each of the 0- 
cycles which are the descendants of the left n-cycle. On 
the other hand, the peak number always maintain to be 
4 for each of the 0-cycles belonging to the descendants of 
the right n-cycle along with the decreasing radius as n 
increases. 



V. JULIA SET AND INFINITE N LIMIT 

The recursion relation given by Eq. (15) happens to be 
a rational map of degree 4. Then, from the work of Julia 
and Fatou, we may conclude that the backward iterations 
defined by Eq. (16) leads towards the Juha set associated 
with the map of Eq. (15). 

Generally, Julia sets can be divided into two classes: 
Some are connected in one piece while the others are just 
a cloud of points. The tendency of the zero density with 
increasing n shown in the above section indicates that 
the Julia set here belongs to the latter. Moreover, by 
observing the distribution patterns and the densities of 
the zeros, we may conclude that a 0-cycle belonging to 
the descendants of the right 71-cycle shrinks to a point 
in the limit of infinite n, and these infinite number of 
points coincide not only with the 0-cycles generated from 
the left n-cycle but also with the zeros obtained from the 
condition tanhry^""^) = ±i for infinite n. Thus, the same 
Julia set arises in DHIM as well as DDIM in the infinite 
limit. In fact, the loci of the zeros of DDIM in the infinite 
limit are identical to that of DHIM. 

The Julia set arising from the distribution of the Fisher 
zeros possesses multifractal structure, and the corre- 
sponding generalized dimensions Dq and singularity spec- 
trum / (a) , obtained by using derivative method and by 
approximating the limiting set of the zero distribution 
with that of n = 8 shown in Fig. 6 [l^ ITEI l , are shown in 
Fig. 7. ^ 

The sinh (2r]) complex plane is also widely used in 
studying the distribution of the zeros for simple square 
Ising model. The results of simple square Ising model im- 
ply that the zeros may lie on the unit circle, |sinh 2?]*^"^ | = 
1, for an n-lattice. To obtain the distribution in the 
sinh (2ry) complex plane, starting with the unit circle we 
then perform n-fold backward iterations of the map. 



sinh 



(2^(")) 



(sinh (2r/("-i)))'^ 
+ (sinh (277("-i)))^ 



2 (sinh (277("-i))) 
1 + (sinh(2?7("-i)))^ 



(40) 



which is an equivalent expression of Eq. (6). This is also 
a rational map of degree 4, and the Julia set associated 
with this map can be approximated with n-fold back- 
ward iterations for a sufficiently large n. The resultant 
distribution of n 8 are shown in Fig. 8. Though the 
distribution pattern is different from that in the complex 
tanh (77) plane shown in Fig. 6, owing to the fact that 
Julia set is an invariant set the global multifractal struc- 
ture characterized by Dq and / {a) is the same as that 
shown in Fig. 7. 

The density near the bulk transition point has a linear 
behavior, and the linear behavior gives the logarithmic 
singularity of the specific heat. By observing the results 
of Eqs. (26), (33), and (38), we may expect that the 
linear behavior for the density near the bulk transition 
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point disappears in the limit of infinite n. This leads to 
the cusp behavior in the specific heat at the critical point 
as the decoration- level goes infinite [23| ■ 

VI. SUMMARY AND DISCUSSION 



In summary, we study the distribution and the density 
of the Fisher zeros for the Ising model defined on square 
lattices with diamond-type bond-decorations in this pa- 
per. By carrying out the exact renormalization map of 
temperature, we can express the free energy of an n- 
lattice as the sum of two parts: One, referred as the local 
part, is mainly the contribution from an n-level decorated 
bond; and the other, referred as the long range part, is 
the contribution from the interactions among n-level dec- 
orated bonds. For the local part, we obtain the zeros as 
a set of scattered points. Along with the increase of the 
number of the scattered points in the distribution of the 
zeros as the decoration-level n increases, all zeros be- 
longing to the lower decoration-levels are also contained 
in the zeros distribution of the higher decoration-level. 
The correponding zeros of the long-range part are the 
union of continuous closed curves which also form sub- 
rings due to the intersections. For the long-range part, 
the distribution of the zeros leads towards the Julia set 
associated with the renormalization map which is a ra- 
tional map of degree 4. This Julia set also serves as the 
limiting set of the zeros obtained from the local part of 
the free energy. The zero representing the critical point 
of ferromagnetic phase transition is one of the zeros of 
the long-range part in the physical region, and the locus 
of the critical point for an arbitrary decoration-level n is 
given. Along with the distribution pattern, we also cal- 



culate the density of the zeros of the long-range part for 
the cases of n = 1 and 2. The evolution of the density 
of zeros in the course of increasing n indicates that the 
Julia set is a cloud of points. Thus, the Julia set has the 
Hausdorff dimension greater than one and possesses the 
multi-fractal structure. 

Comparing with the continuous closed curves appear- 
ing in the zeros distribution of the simple square Ising 
model, we have more complicated distribution patterns 
for the bond-decorated Ising model. But, the patterns re- 
main to be continuous closed curves in both tanh (77) and 
sinh (2r]) complex planes for a finite decoration-level n. 
This may be due to the fact that the interactions among 
the Ising spins are effectively isotropic after the renor- 
malization map to the corresponding 0-lattice, although 
the coordination numbers of lattice sites are highly in- 
homogeneous for an n-lattice with large n. It has been 
demonstrated on classic lattices that only when the cou- 
plings among the nearest neighbors change from isotropic 
to anisotropic, the distribution of the Fisher zeros may 
change from continuous curves to an area in the plane 
20] . This leads to the conclusion that in the limit of 
infinite decoration-levels the system has completely dif- 
ferent properties: The distribution of the zeros has multi- 
fractal structures, and the nature of phase transition of 
the system is different from that of finite decoration-levels 
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fig.l A decorated bond with the decoration level (a) 
n = 0, (b) n = 1, and (c) n = 2. Note that the Ising 
spins with the Latin subscripts are referred as inner spins. 

fig. 2 (a) The distribution of the Fisher zeros of an n- 
lattice in the complex tanh plane as the left (gray) 

and the right (bold-faced) n-cycle, and (b) the density 
(left vertical scale) and the raidus (right vertical scale) 
of the Fisher zeros in the left n-cycle. 

fig. 3 (a) The four (n — l)-cycles for the zero distribu- 
tion of an n-lattice in the complex tanh (77^"^^^) plane. 
The cycles displayed by the bold-faced curves are the de- 
scendants of the right n-circle and the gray curves are 
those from the left n-circle. The points indicated by 
crosses are taken as the centers of the (n — l)-cyclcs, and 
they arc the prcimages of the renormalization map for 
the centers of two n-cycles. (b) The density (left vertical 
scale) and the radius (right vertical scale) of the Fisher 
zeros in the (n — l)-cyclc with the center at (1,0). 

fig. 4 (a) The 12 (n — 2)-cycles for the zero distribu- 
tion of an n-lattice in the complex tanh (ry^"~^^) plane. 
The points indicated by crosses are taken as the centers 
of the (n — 2)-cycles, and they are the preimages of the 



renormalization map for the centers of the (n — l)-cycles. 
The 12 (n — 2)-cycles can be divided into three classes, 
/, II, and ///, according to the decreasing order in the 
magitude of the circumference of a circle. The density 
(left vertical scale) and the radius (right vertical scale) of 
(b) the (n — 2)-cycle, which is the one marked / in (a), 
with the center coordinate (1,0), (c) the (n — 2)-cycle, 
which is the one marked II in (a), with the center co- 
ordinate vl+Vle*'^/*, and (d) the (n — 2)-cycle, which 
is the one marked III in (a), with the center coordina 

fig. 5 The zero distribution of an n-lattice in the com- 
plex tanh (77) plane obtained from the condition of Eq. 
(12) for n = 4. 

fig. 6 The zero distribution of an n-lattice in the com- 
plex tanh (77) plane obtained from the condition of Eq. 
(12) for n = 8. 

fig. 7 The (a) generalized dimension Dq and (b) sin- 
gularity spectrum / (a) of the Julia set associated with 
the renormalization map. 

fig. 8 The zero distribution of an n-lattice in the com- 
plex sinh (2?7) plane obtained from the condition of Eq. 
(12) for n = 8. 
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